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The concept of Schramm-Loewner evolution provides a unified description of domain boundaries of many 
lattice spin systems in two dimensions, possibly even including systems with quenched disorder. Here, we 
study domain walls in the random-field Ising model. Although, in two dimensions, this system does not show 
an ordering transition to a ferromagnetic state, in the presence of a uniform external field spin domains percolate 
beyond a critical field strength. Using exact ground state calculations for very large systems, we examine ground 
state domain walls near this percolation transition finding strong evidence that they are conformally invariant 
and satisfy the domain Markov property, implying compatibility with Schramm-Loewner evolution (SLE re ) with 
parameter n — 6. These results might pave the way for new field-theoretic treatments of systems with quenched 
disorder. 



In the past decades, analytic techniques such as conformal 
field theory (CFT) and Coulomb gas methods have led to a 
rather comprehensive understanding of critical phenomena in 
two dimensions (2D). In particular, CFT allows for a com- 
plete classification of 2D critical points, the exact determi- 
nation of critical exponents and, in some cases, even scaling 
amplitudes [ 1 1. This success is tied to the fact that the confor- 
mal group is infinite-dimensional, however, which is true only 
in 2D, and few of the results generalize to higher dimensions 
0. Another difficulty for this approach arises for the impor- 
tant class of systems with quenched disorder, such as diluted 
magnets, random-field systems and spin glasses [3 |, since the 
non-unitary CFTs that are believed to describe systems with 
quenched disorder are poorly understood [4]. 



While some geometrical aspects of critical phenomena had 
been previously worked out using concepts from the Coulomb 
gas Q and two-dimensional quantum gravity (6), a break- 
through was achieved with the description of domain bound- 
aries in terms of random curves in the plane in a frame- 
work dubbed Schramm-Loewner evolution (SLE) Q. In 
SLE, stochastic curves in the plane are constructed from one- 
dimensional Brownian motion, thus classifying a statistical 
ensemble of curves with only one parameter, the diffusion 
constant k. Characteristic interfaces in many physical sys- 
tems have been shown (in some cases rigorously) to satisfy 
SLE K . These include percolation (k — 6), self avoiding 
walks (k = 4/3), as well as spin cluster boundaries (k — 3) 
and Fortuin-Kasteleyn cluster boundaries (k — 16/3) in the 
Ising model. In recent years, close connections between SLE 
and CFT, including links between probabilistic properties of 
SLE curves and scaling operators in CFT, or between the cen- 
tral charge c of the CFT and the diffusion constant k have 
been established [7]. A number of numerical studies have 
found interfaces in disordered systems to be (partially) con- 
sistent with SLE, in particular the 2D Ising spin glass El ED, 
the Potts model on dynamical triangulations fit)! , the ran- 
dom bond Potts model ifTTl . and the disordered solid-on-solid 
model [ 12 1. Such findings and the close link between SLE and 
CFT spur the hope of a more complete understanding of sys- 
tems with quenched disorder from a field-theoretic perspec- 
tive. 



I. RANDOM-FIELD ISING MODEL 

Here, we study domain walls in the random-field Ising 
model (RFIM) with Hamiltonian fl21 

"K = - J} j SjSj - J] hjSj, (1) 
{id) 1 

where the spins Sj = ±1 are located on the sites of a square 
lattice and interact ferromagnetically with nearest neighbors. 
The local fields hi are quenched random variables which, for 
the time being, we take as drawn from a Gaussian distribution 
with mean H and standard deviation A. Since only the ratio 
J I A is relevant, we take J = 1. Random field models have a 
large number of experimental realizations which are of tech- 
nological importance such as superfluid helium, liquid crys- 
tals in silica aerogels, Bragg glasses in high-T c superconduc- 
tors, amorphous solids| 14- 1 and ferroelectric materials [ 15 1. It 
was shown by Imry and Ma [ 16] that random fields destabi- 
lize the ferromagnetic order in dimensions d < 2. For the 
case of 2D, it was argued by Binder IfTTl that ferromagnetic 
order occurs only up to a break-up length scale Lb ~ e A / A , 
which grows with decreasing disorder A, and that the system 
remains paramagnetic at scales L > Lf,. Later, Aizenman 
and Wehr 1 1 8 1 proved that for d < 2 the system indeed has 
a unique Gibbs state, precluding the existence of an ordering 
transition. On the contrary, for d > 3, diverges at the ther- 
modynamic transition point, below which the system is ferro- 
magnetic 1T31 . For non-zero average fields H, on the other 
hand, even in 2D the size of spin clusters diverges at a critical 
value H c = H C (A) |fi"9H22"l . However, the weight of these 
clusters is sub-extensive, such that the free energy remains 
analytic and no thermodynamic phase transition occurs. This 
phenomenon bears many similarities to the Kertesz line in fer- 
romagnets in the absence of disorder ll23ll . It is this transition 
at non-zero H that we study in this Letter. 

To investigate the properties of domain walls in the RFIM 
we numerically compute exact ground states of samples of 
random-field realizations. Ground states can be found in poly- 
nomial time via a mapping to a minimum-cut/maximum-flow 
problem 11241 . We employ a fast algorithm based on the idea 
of "augmenting paths" 11251 . which allows us to find ground 
states of systems of 10 7 spins in about 6 s, such that the maxi- 
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mum system sizes exceed those of previous studies |[T9lj2Tl by 
about an order of magnitude. We use a variety of domain ge- 
ometries, partially with fixed spins to enforce the occurrence 
of domain walls. The calculations reported here were carried 
out at either A = 2.5 and H = 0.01362 = H C (A) ± 0.00007 
or at A = 1.7 and H = 5.08 x 10~ 4 = H c (A)±0.07x 10~ 4 . 
For both cases the breakup length scale is only a few lattice 
spacings, much less than the system sizes we look at. Both 
the breakup length scale and the critical external field were 
determined using recipes laid out in Ref. [20]. 



II. SCHRAMM LOEWNER EVOLUTION 

In the framework of Loewner evolution one imagines a ran- 
dom curve jt in the plane as being continuously grown in time 
t in a random process. Instead of studying this process di- 
rectly, one considers the evolution of a family g t : H \ -f t —> H 
of conformal maps that take the complement of -f t in the upper 
half plane H to H. Under this map the curve j t , which lies on 
the boundary of H \ j t , is taken to the boundary of H, i.e. to 
the real line. Assuming standard normalization and boundary 
conditions, it turns out that g t is completely determined by the 
one-dimensional function £ t , which corresponds to the image 
under g t of the tip of the growing curve on the real line, via 
the differential equation 

dgt(z) = 2 

It was shown by Schramm [26] that if the ensemble of curves 
7 f is conformally invariant and satisfies the domain Markov 
property (to be discussed below), the one-dimensional random 
process described by £ t must be Brownian motion with zero 
mean and variance nt. For such SLE curves, many stochas- 
tic properties can be calculated rigorously including, for in- 
stance, the fractal dimension df = 1 + k/8 or the probabil- 
ity Plp(2;, v) that the curve -f t passes to the left of the point 
(x,y). The latter was proven by Schramm [27 1 for curves 
starting at the origin of the upper half plane H to be 

(3) 

where 2-F1 is Gauss' hypergeometric function. 



III. DOMAIN MARKOV PROPERTY 

The domain Markov property (DMP) formalizes the notion 
that a growing path of the type described above is agnostic 
about its past. Let Po("f a b) be the probability measure of 
curves j a b in a domain D running between points a and b on 
the boundary of D, and let c be a point in the interior of D. 
Then, the DMP states that 

Po(lcb\lac) = Po\y ac (lcb), (4) 




FIG. 1. (Color online) The main panel shows the spatial distribution 
of \APomp{x, y)\ for a system of 4000 spins in a circular domain 
and a vertical cut 7„ c of length 6 starting at the bottom of the do- 
main. Lighter colors correspond to larger deviations from the DMP. 
The right inset panel shows the decay of j APdmp(35, y)\ with dis- 
tance in the positive vertical direction Ay (in lattice spacings) from 
the point of maximum deviation. Three different system sizes and 
correspondingly scaled cut lengths 7 ac are shown. The left inset 
panel shows the decay in the negative vertical direction (i.e. along 
the cut) vs. Ay scaled by the cut length. The tests were performed at 
A = 2.5 and H = 1.362 x 10" 2 « H C (A). 

i.e., the probability of j co is independent of whether 7 ac is 
preconditioned in domain D, or whether 7 ac is excluded from 
the domain itself. While it is not debated that the DMP holds 
for domain boundaries in pure lattice systems Q, even off 
criticality, it has been argued that it is likely not to survive the 
average over quenched disorder [9|. For the RFIM, we have 
checked the DMP numerically using the left passage proba- 
bilities Plp(x, y) instead of calculating the probability of all 
possible curve segments -f c i,. For the l.h.s. of Eq. this 
amounts to picking out those configurations of the random 
fields that yield an interface along 7 ac , while for the r.h.s., 
the interface is asserted to run along j ac , for instance by fix- 
ing the corresponding spins with large magnetic fields. If the 
DMP holds, then 

APdmpO?, y) =y^ Pa(jcb\7ac)PLp{x, y; icb) 

^ (5) 

- 2^ Pa\ 7a Slcb)Php(x, y; j c b) 

leb 

will be identically zero. We have studied system sizes of 1000, 
4000, and 16 000 spins with proportionally scaled cuts 7 ac 
of length 3, 6, and 12, respectively. For the largest system, 
we looked at 3 x 10 s ground state configurations, of which 
only about 2800 satisfied the conditioning along -f ac . As is 
clearly seen from |APdmp(^, y) \ shown in the main panel of 
Fig.[T] the DMP does not survive the disorder average exactly. 
The deviations are maximal around the tip of j ac and fall off 
rapidly with distance from j ac . As shown in the right inset 
panel the decay of APqmp with the vertical offset Ay from 
the tip of -f ac is nearly independent of the system size and the 
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length of 7 ac . In contrast, as shown in the left inset, the decay 
rate for Ay < 0, i.e., along the cut, is proportional to the cut 
length. Perpendicular to the cut, the decay rate (not shown) 
is again independent of system size. Hence, the intrusion of 
deviations into the interior of the domain extends to only a few 
lattice spacings and is largely independent of system size and 
cut length, such that the DMP will be recovered in the scaling 
limit. We find that the agreement with the DMP in the scaling 
limit also holds off the critical percolation line. 

IV. LEFT PASSAGE PROBABILITY 

We examined the agreement of the RFIM interfaces with 
the SLE expectations for the left passage probabilities. As the 
rigorous result of Eq. ^ is valid on the upper half plane H, 
we performed our ground state calculations on lattices embed- 
ded in domains D which have simple, closed-form conformal 
maps w(z) to H [12]. By fixing the boundary spins through 
the respective random fields, the interface was forced to run 
between points a and b in D which are mapped to the origin 
and infinity in H, respectively. Numerical checks were per- 
formed for the unit circle with the interface between —i and 
i, which is mapped to H by w(z) = Jyrf , as well as the unit 
square with interfaces defined from to l+i, which is mapped 
to H by w(z) = — p(l +i — z; 1, i) where p(z; wi,w 2 ) is the 
Weierstrass p-function. Looking at the left passage probability 
for multiple domains additionally acts as a check of conformal 
invariance. For a quantitative comparison we considered the 
mean square deviation of the computed left passage probabil- 
ity Plp from the exact result P£ F of Eq. 

E(k) = ({P LP (x,y) - P£ P (x,y)] 2 y^ , (6) 

where (-)d denotes a spatial average over D, excluding the 
vicinity of the fixed boundary spins. This quantity is shown 
in Fig. [2] as a function of k for the circle and square geome- 
tries. In both cases E(k) is minimal for k within 0.05 of the 
value k = 6. The spatial dependence of the deviation from 
P£p 6 (x, y) is also shown in Fig. [2] and, for a horizontal cut 
through the domain, in Fig. [3] There appears to be no system- 
atic deviation. 



V. CROSSING PROBABILITY 

From the presented results it is evident that spin domain 
interfaces in the 2D RFIM satisfy the domain Markov prop- 
erty and are conformally invariant at the percolation threshold 
H = H C (A) in the scaling limit, and thus are described by 
SLE. Furthermore, the parameter k appears to be consistent 
with k — 6 to high precision. To corroborate these findings, 
we tested our results for compatibility with the exact formu- 
las for the crossing probabilities of percolation clusters, an- 
other system with k = 6 l28l . These give the probability 
of finding a cluster of a given species (say up spins) which 
touches two non-adjacent segments of the boundary of a do- 
main. In particular, the probability ir r of a domain touching 
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FIG. 2. Magnitude of deviation of the left passage probability for 
spin domain interfaces from the exact result of Eq. {3} for n = 6 in 
circular and square domains. The interfaces are constrained to run 
between points a and b as shown. The lower panel displays the spa- 
tially averaged deviation E(k) as a function of the diffusion constant 
k, showing a clear minimum close to n — 6. Calculations were per- 
formed for 10 000 disorder realizations on systems of 6 x 10 6 spins 
at A = 2.5 and H = 1.362 x 10~ 2 « H C (A). 
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FIG. 3. Deviation from the exact left passage probability along a hor- 
izontal line crossing the center of the circle geometry that is shown 
in Fig. [2] 



both the top and bottom of a rectangle of aspect ratio r is 
known to depend only on r l28l . Similarly, in a domain de- 
fined by an equilateral triangle, the probability of a cluster 
crossing from a fraction x of one boundary edge to the oppo- 
site edge is tt x = x ll29l (see Fig.[4]for a schematic representa- 
tion). Figure |4] shows these exact percolation results together 
with numerical simulation data for the RFIM at A = 1.7 and 
H = 5.08 x 10~ 4 rj H C (A = 1.7) for a system of 6 x 10 6 
spins for both the rectangle and triangle geometries. We find 
very good agreement with the percolation results, cf. the lower 
panel of Fig. [4] We also show results calculated with an exter- 
nal field H = 4.71 x 10~ 4 w 0.93 H c . Systematic deviations 
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FIG. 4. Upper panel: Crossing probabilities for spin clusters in 
the 2D RFIM in rectangular (lower scale) and triangular (upper 
scale) domains. Calculations were performed at A = 1.7 and 
H = 5.08 x 10" 4 w H c (A) for20-30xl0 3 disorder realizations 
for the rectangular domain and 7xl0 3 realizations for the triangu- 
lar domain. The solid lines indicate the exact results for percolation 
clusters (SLE K= 6) derived in Refs. [28 ,29]. Lower panel: Deviations 
of the numerical results from the exact expression for percolation for 
the rectangular domain. Deviations in the triangular domain are sim- 
ilar. Shown with dashed error bars are data calculated at A = 1.7, 
but at an external field H = 4.71 x 10~ 4 . Some systematic devia- 
tions from SLE expectations are already seen for this slight detuning 
from criticality. 



from SLE can be clearly seen, even for this slight detuning 
from criticality, indicating the high sensitivity of our tests. 



VI. FRACTAL DIMENSION 

One of the rigorous results for curves described by SLE K is 
their fractal dimension, which is given, for k < 8, by df — 
1 + k/8. We numerically determined the fractal dimension 
for a range of different geometries and display the results in 
Fig. [5] We find that corrections to scaling for the interface 
length Li are well described by the form 



aL d f(l + b/L). 



(7) 



We take L as the square root of the number of spins for the 
case of non-rectangular domains. For the circle geometry we 
find df — 1.7506(9), while for the square geometry we arrive 
at d f = 1.7514(14), using A = 2.5 and H = 1.362 x 1CT 2 
H C (A) in both cases. These results are perfectly compatible 
with df — 7/4 = 1.75 expected for SLE curves with k = 6. 
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FIG. 5. Scaling of interface lengths with system size. The best 
fit values — df = 1.7506(9) for the circle geometry and df = 
1.7514(14) for the square geometry — are in excellent agreement 
with the value expected for SLE6, df = 7/4. The calculations were 
carried out at A = 2.5 and H = 1.362 x 10" 2 « H C (A). 



VII. BROWNIAN MOTION 

As the most direct test for SLE, we studied the one- 
dimensional stochastic process (or driving function) £ ( gen- 
erated by the Loewner map g t according to Eq. |2]) as applied 
to domain walls in the RFIM. For a lattice system, the family 
of maps g t is realized as a discrete series of maps gi iteratively 
removing a small section from the beginning of the curve. For 
this purpose, g L is approximated using a vertical slit map 11301 



9i(z) = iy/-(z- £i) 2 -4AU + & 



(8) 



Here, £, and Ati are determined through £j = Xi t i-i and 
At, = yf j_i/4, where x^j-i and are the coordinates 

of the i'th segment of the curve after undergoing the i — 1 suc- 
cessive maps g^_i o . . . o g 1 . The parameter & is the value of 
the driving function £ t sampled at time ti = J2j<i &tj. The 
complex square root in Eq. (|8]l is calculated, as usual, with 
the branch cut along the negative real axis. We studied the 
statistics of 10 000 interfaces generated in a half disc, opti- 
mally mimicking the full space H. The interface is initiated at 
the origin by two fixed spins and is considered ended when it 
touches the curved boundary. We used systems of 6 million 
spins at A = 2.5 and H = 1.362 x 10~ 2 « H C (A). We find 
that the variance of the driving function calculated from the in- 
terfaces is k = {{^ - (6)) 2 ) A = 6.086(87), and the normal- 
ized mean is £ = (£t)/v&i = 0.017(10), perfectly compati- 
ble with SLE K= 6. Using a Kolmogorov-Smirnov test 113011311 . 
we further checked that the £j are normally distributed and 
find a p-value of p = 0.17, indicating consistency with a nor- 
mal distribution. To check for the statistical independence of 
the increments of £ t , we divided £ t into n increments evenly 
spaced in time and checked whether the signs of these incre- 
ments follow a x 2 distribution with 2" — 1 degrees of freedom 
ll30ll . For n — 10, we find a p- value of p — 0.18 indicating 
consistency with the assumption of statistical independence. 
As a further check of conformal invariance, we performed 
the same tests on interfaces originating in the two domains 
of Fig. [2] and found similar agreement with SLE. 
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FIG. 6. Crossing probabilities for the binary distribution of fields 
as a function of the fraction p of up fields. The upper panel shows 
the dependence of these curves on disorder strength (for system size 
L 2 — 10 6 ). The bottom panel shows the crossing of these curves 
for different system sizes at the critical point p c (at disorder strength 
A = 2.1). 



VIII. BINARY FIELDS 
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FIG. 7. Left passage probabilities for the binary distribution of 
fields. See Fig. [2] for a detailed description of the setup. The cal- 
culations were carried out at A = 2.1 and p = 0.522 « p c (A) for 
systems of 6 million spins. Averages were performed over 10 ran- 
domly chosen ground states for each of 6000 disorder configurations. 



demonstrating that different sizes cross at p c . We also point 
out that this crossing of curves happens when the crossing 
probability is 0.5 as expected. Interestingly, as can be seen 
in in the upper panel of Fig. [6j for A < 2 this transition 
appears to occur at the constant value p c (A < 2) = 1/2, 
while p c (A > 2) > 1/2. We test agreement with SLE 
predictions by looking at the fractal dimension at p c , finding 
df = 1.746(2), and by looking at the left passage probability, 
the results of which are shown in Fig. [7] Both properties are 
consistent with SLE for k = 6. 



Finally, we also considered a binary (Bernoulli) field dis- 
tribution, where each local random field hi takes on the value 
A with probability p and —A with probability 1 — p. Be- 
cause of the discrete nature of the distribution, this system has 
a massive ground-state degeneracy, and behaves rather differ- 
ently from the Gaussian RFIM, at least at zero temperature. 
Although polynomial algorithms for enumerating all ground 
states are known |32], handling all ground states becomes im- 
practical for larger system sizes. Here, instead, we sample 
from the ground-state manifold by adding a tiny noise term 
(normally distributed with strength 8) to the Bernoulli field 
distribution. For sufficiently small 8, the resulting state is 
also a ground state of the noise-less system, and, importantly, 
is selected without bias from among the degenerate ground 
states. We find that there exists a geometric transition where 
spin clusters diverge at p = p c (A), in analogy with H C (A) 
for the Gaussian case. This is illustrated in Fig. [6] where 
we have plotted, for a number of different disorder strengths, 
the probability of finding a spin-up cluster that touches both 
the top and bottom boundaries of a square geometry. Shown 
in the bottom panel is the size dependence of these curves, 



IX. CONCLUSIONS 

We have studied the properties of spin cluster interfaces in 
the ground state of the Gaussian random field Ising model at 
values of the external field strength H where the size of the 
clusters diverges. For this T = system with quenched dis- 
order, the domain Markov property was shown to be satisfied 
in the scaling limit. Together with the conformal invariance 
of the interfaces deduced from Schramm's formula and the 
crossing probabilities, it is shown clearly that the spin domain 
interfaces satisfy SLE K=6 , corresponding to pure percolation. 
The fractal dimension is in perfect agreement with these ob- 
servations, contrary to the case of the solid-on-solid model 
studied in Ref. iTPZl . where k ~ 4 was found from Schramm's 
formula, but dt ~ 1.25 ^ 1 + re/8. Studying the SLE map 
directly, we have shown that the driving function describes 
Brownian motion. The consistency with SLE carries over to 
the case of binary random fields, where degeneracies occur. 
This is in contrast to the observations for the spin glass model, 
where domain walls appear to be only described by SLE for 
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continuous disorder distributions |00[9l[33|. The 2D RFIM 
thus seems to provide a paradigmatic example where SLE is 
realized in all known aspects in a system with quenched dis- 
order, nourishing the hope for a more systematic treatment of 
systems with quenched disorder in field theory. 
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